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The giant magneto-resistance (GMR) of ferromagnetic bilayers with a supercon- 
ducting contact (F1/F2/S) is calculated in ballistic and diffusive regimes. As in 
spin-valve, it is assumed that the magnetization in the two ferromagnetic layers Fl 
and F2 can be changed from parallel to antiparallel. It is shown that the GMR 
defined as the change of conductance between the two magnetic configurations is an 
oscillatory function of the thickness of F2 layer and tends to an asymptotic positive 
value at large thickness. This is due to the formation of quantum well states in 
F2 induced by Andreev reflection at the F2/S interface and reflection at F1/F2 in- 
terface in antiparallel configuration. In the diffusive regime, if only spin-dependent 
scattering rates in the magnetic layers are considered (no difference in Fermi wave- 
vectors between spin up and down electrons) then the GMR is supressed due to the 
mixing of spin up and down electron-hole channels by Andreev reflection. 

PACS: 74.80.Dm, 75.70.Pa, 75.70.Cm 

The mechanisms of Giant Magnetoresistance (GMR) in magnetic multilayers and sand- 
wiches are well established. They are related to the influence of the spin assymetry on 
the conductivity. Several microscopic factors may play a role: the electronic band struc- 
ture (e.g. exchange splitting) Q has a direct influence whereas the spin assymetry in the 
bulk or interface electron scattering rates has an indirect influence [fj. Experiments on 
GMR are carried out in two geometries: current in plane of the structure (CIP) and 
current perpendicular to plane (CPP) [|J. It was observed that the CPP-GMR amplitude 
is always larger than the CIP one [[| . It should be mentioned that in CPP measurements, 
superconducting leads are most often used as contacts on the spin- valve structure jJJ. 
Recently, in Ref. || the CPP-GMR of Co/Cu layer was numerically calculated using a 
realistic band structure for Co and Cu and in the presence of one superconducting contact. 
Unexpectedly, it was found that in this case the GMR is completely supressed due to An- 
dreev reflection on the ferromagnet /superconductor (F/S) interface. In order to resolve 
this contradiction between theory and experiment, we developed an analytical theory of 
CPP-GMR for a spin valve sandwich of the type F/F/S, where F's are ferromagnetic lay- 
ers, the magnetizations of which can be oriented parallel or antiparallel to each other, S 
is a superconducting contact. A simple two band (spin up and down) free electron model 
is adopted for this calculation. 

For calculating the conductance of the considered system, we used the generalized 
Fisher-Lee formulae in spinor form ||: 



a a ^\z } z') 



where V 2 = ~ * s ^ e antisymmetric gradient operator and 

Af(z,z') = (i/2)(Gf(z,z r ) - GP a *(z',z)). Here G aa is the conventional Green function 
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and G af3 (a 7^ (3) is the anomalous Green function, antisymmetric in spin indices (oc,/3). 
The summation is performed over available channels, k = (k x ,k v ) is the component of 
electron momentum in the XY-plane of the layers and z is the coordinate perpendicular 
to the XY-plane. 

The Green functions satisfy the system of Gor'kov equations 0: 

(«w - ~ e F ) + e ex (z))G^(z, z') + A(z)G^(z, z') = 5{z - z') 

(ico + - e F ) + e ex (z))G^(z, z') + A*(z)G^(z, z') = 0, (2) 

where the superconductor gap A(z) is considered as constant inside the superconductor 
and zero in the ferromagnetic layers. On the contrary, the exchange splitting parameter 

e ex (z) = ^^(k^ F — k^ F ) is zero in superconductor and constant in the ferromagnetic layer 

(kp represent Fermi momenta for spin up (down) electrons ). The system of equations 
(2) can be solved exactly in the clean limit: e.g. if the thicknesses of the ferromagnetic 
layers are much smaller than both the elastic mean free path and the magnetic length 

4=^-, where D is the diffusion constant. The same assumption was made in the first 
part of reference 0. 

The expressions of the conductances in parallel (P) and antiparallel (AP) alignments 
of the magnetization in the adjacent ferromagnetic layers can be written in the following 
form 

e^£(i-« 2 ), (3) 

G"=^E ' 1 -^' 1 - r2 > , (4, 

TlTl ^ l+2r 2 R 2 +r 4 -2Rr(cos2c^a-cos2c 1 a)(l+r 2 )-2r 2 (R 2 cos 2(cT+ c J-)a-cos 2(cT -c L )a) 
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where r = ^ — represents the effective spin polarization, R = -^-^ 1, c^' = 

c' + & ' c'c 1 + c 2 

Re\J (kp') — k 2 , and c 2 = ReJ (k s F ) 2 — k 2 . a is the thickness of the intermediate fer- 
romagnetic layer. k s f is the Fermi wave vector in the superconducting layer. The upper 
limit of the sum over n is equal to the minimum value of k F or kp. 

The physical meaning of the obtained expressions is rather clear. In the P-configuration, 
the conductance of the system decreases compared to its value in the absence of super- 
conducting contact, due to Andreev reflection. This conclusion coincides with the result 
of the numerical calculations of 0. If r = 0,( e.g. the paramagnetic metal in contact with 
superconductor), expression (3) coincides with known results (See eq. (127) in Ref. ||). 
Expression (4) contains two factors (1 — R 2 ) and (1 — r 2 ). The first one is due to Andreev 
reflection on F/S interface and the second one is the usual reflection of electrons at F^/F^ 
interface. So, if only these two factors are taken into account, considering the denomina- 
tor in (4) equal to unity, a finite GMR amplitude is still obtained. Let's consider then 
the effect of denominator in expression (4). It describes the multiple reflections of an 
electron which moves inside the ferromagnetic layer adjacent to the superconductor, as in 
a Fabri-Perro interferometer. These multiple reflections are responsible for the formation 
of quantum well states within the layer. As a result, the conductance G AP is an oscillatory 
function of the arguments k F a, (k F ± kp)a, but it never diverges nor becomes negative. 
A similar behaviour of conductance was predicted in |§, for a structure composed of a 
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ferromagnetic layer sandwiched between a superconducting contact on one side and a thin 
oxide barrier on the other side. 

Now we come to the question whether G p is always larger than G AP , or if for some 
values of parameters, G AP can be equal or even larger then G p . For very small polarization 
r « 1, from expressions (3) and (4), it is easy to obtain the following approximate 
expression for the GMR: 

G M«=^ = 2 ? ^ = 2 ? (^) 2 (^|) ! (5) 

This expression is definitely positive. Without any superconducting contact, the GMR 
would be given by GMR = 2J2 r2 - R is interesting to note that expression (5) coincides 

K 



with the expression of the MR in a spin- valve tunnel junction [10], after substitution of 
C2 by the modulus of the imaginary electron momentum inside the barrier. The physics 
of both phenomena is similar: in both cases the electrons undergo reflections on the 
interface: F/I (I-insulator) or F/S. These spin-dependent reflections change the spin- 
dependent density of states in the ferromagnet near the interface and, correlatively change 
the polarization of the current. 

For larger r, the conductances G p and G AP are plotted in fig. 1 versus the square of the 
effective polarization r 2 . The following parameters were used: k F = 1A~\ kp = 1.3A" 1 
and a = 5co, (cq = 4.06A is the lattice parameter of Co for hep structure), a = 5co 
corresponds to 10 atomic monolayers, kp was varied from 1 to 0, so that correspondingly 
r 2 was changing from to 1. As can be seen in fig. 1, the conductances for both magnetic 
configurations decrease as r 2 increases. G AP exibits also some weak oscillations as a 
function of r 2 , but it remains smaller than G p for almost the whole range of r 2 . For 
r 2 > 0.9, the GMR saturates at a value ~ 220% (see fig. 3) but it does not diverge if 

2r 2 

r — > 1 as it is the case for a spin-valve without superconducting contact (GMR = — - — 2 
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Figure 2: The difference of the conductivities Aa = op — a^p for the parallel and an- 
tiparallel configurations as a function of the thickness of the ferromagnetic layer a. The 
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parameters are: k s F = 1.3 A , k F — 1.0 A 
correspond to integer numbers of atomic layers of Co. 
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in this case). 

Now let us look at the effect of averaging over a distribution of magnetic layer thickness 
variations since it could arise that these fluctuations may suppress the GMR. In fig. 2, the 
difference G p — G AP is plotted versus the thikness a for a given value of r 2 = 0.16. Fig. 2 

qP _ qAP 



shows that this difference, and consequently the GMR 
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oscillates around a 



non-zero positive value and tends to the asymptotic limit 44% for a > 100A. Of course, 
the thickness of the layer can change only by steps equal to the lattice parameter Co/2. 
In fig. 2 the possible values of a have been marked considering Co = 4.06A for Co. It is 
interesting to note that the situation is similar to the case of a spin-valve tunnel junction 
with paramagnetic metal layer inserted between one ferromagnetic electrode and an in- 
sulating barrier. ||11|| . In this case, it was shown that the paramagnetic layer (for instance 
Cu inserted between Co and AI2O3) can constitute a spin-dependent quantum well. Oscil- 
lations in tunnel magnetoresistance (TMR) were predicted for such system as a function 
of the paramagnetic layer thickness with a period given by the Fermi-wave length in this 
layer. However, a crucial difference between this case and the present one is that here, 
the GMR oscillates around a finite positive value whereas in a tunnel junction, the GMR 
oscillates around zero. Consequently, for tunnel junctions, averaging over a distribution 
of the paramagnetic layer thickness caused by roughness, and/or increasing the param- 
agnetic layer thickness leads to a strong decay in TMR amplitude. In contrast, in the 
present case, averaging over a distribution of thickness and/or increasing the thickness of 
the ferromagnetic layers leads to a non-zero GMR amplitude which depends on the values 
r 2 and R 2 . 

This situation is illustrated in fig.3, where the dependence of the GMR on the effective 
spin polarization r 2 is plotted for two different cases, i.e. with thickness a equal to 10 and 
500 monolayers of Co. In fig.4, the same dependence is presented but for a structure in 
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Figure 3: The GMR as a function of the square of the effective spin polarization r 2 for 
a = 20.3 A - dashed line and a = 1015 A - solid line. Parameters are: kp = 1.3 A -1 , 
k F = 1.0 A" 1 . 

which the layer thikness a is supposed to take random values equally distributed between 

9, 10 or 11 monolayers of Co. We have also calculated the a-dependence of conductivity 

for other values of r 2 . These dependences are similar to the one shown in fig. 2. The 

resonances of G AP exibit rather sharp peaks at a = if 1 — r 2 <C 1, i.e k F — > 0. In 

kp 

this case, the system becomes a real Fabri-Perro interferometer for electrons. 

Now let's consider a different model of the CPP-GMR in spin-valve structures, devel- 



oped in details in [12] and often used for the interpretation of CPP-GMR experiments. 
In this model, it is considered that charge carriers in ferromagnetic metals are s-like 
electrons with negligibly small exchange splitting but with different elastic mean free 
paths for spin up and spin down electrons. The scattering of s-electrons is considered as 
mainly due to s — d scattering, so that the inverse life times of up and down electrons 
are j^L.. = 7 2 cz/] < -^, where 7 is the s — d scattering potential, c is the concentration of 

impurities and v\ 7^ v\ ^> v s are the densities of states states of j and | spin <i-electrons, 
u 8 is the s-electron density of states. Furthermore, we assume that L <^ l s f, where L is 
the thickness of the ferromagnetic layers, l s f is the spin-flip mean free path of s-electrons. 
Now in Gor'kov equations (2), iuo has to be changed to iuj' = iu — ^cG^^ and A 
to A + 7 2 cG^. It is easy to show that if there are no ci-states in superconductor, then 

Gdd, ~ ^ ~ Pd in the ferromagnetic layer, due to simultaneous influence of 

proximity effect on and s — d hybridisation. Consequently, in the following, G dd is 
neglected. It is easy to solve Gor'kov equations and calculate the conductance by using 
expression (1). We have to add to (1) vertex correction but, as it was shown in [JT3 



inclusion of vertex correction is equivalent to a special choice of effective internal electro- 
chemical fields in such a manner that the condition ^ = is satisfied (j is the current 
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Figure 4: The GMR versus the square of the effective spin polarization r 2 for the sandwich 
with the random (9,10 or 11) numbers of Co monolayers. Parameters are: k s F = 1.3 A -1 , 
4 = 1.0 A" 1 . 



in the z-direction). Following this procedure, we found that resistances of the considered 
sandwich with superconducting contacts for parallel R p and antiparallel R AP alignment 
of magnetizations in adjacent F-layers are equal: 

K = R^ P =(a + b)^ + p')/A (6) 

where a and b are the thicknesses of the ferromagnetic layers and p^-' are the resistivities 
for | (|) spin s-electrons. On the other hand, if the S-contact is in a normal state, we 
have 

P _ (a + b)p^p l 

N ~ t i I 
P l + p l 

AP _ (ap T +frp i )(fep T + ap l ) 

(a + 6)(pT+pl) [n 

Therefore within the assumption that the GMR originates from spin-dependent scattering 
rates in the magnetic materials, we find that there is no GMR effect in presence of 
superconducting contact. This conclusion coincides with the results obtained in |L4 |. 



The absence of GMR in this case can be qualitatively understood as follows. In a 
ferromagnetic metal, currents for up and down spin electrons are not equal. However, 
in a BCS superconductor, the current is driven by spin-less Cooper pairs, so that up 
and down spin currents are equivalent. To maintain this equivalence, electrons undergo 
Andreev reflection at the ferromagnet /superconductor interface and spin accumulation 
appears at this interface. Due to this accumulation, a jump Av of chemical potentials 
(of different signs for up and down spin electrons) arises. The values of these jumps are 

Av^ = v^j — P-r = — Av ^ for P configuration and Av^ = —Av- 1 - = v^- — ttt~j — ^tt for 
p T + p l (a + 6)(p T +p^) 

AP configuration, where v is the voltage drop across the total structure. In particular, 
Afj = Ai)[ = for a = b. It can be shown that this additional drop of voltage in parallel 
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configuration exactly equalizes the resistances for P and AP configurations, so that the 
GMR is suppressed. 

In conclusion, contrast to ||, we have shown that in general, due to the exchange 
splitting of electron bands in ferromagnetic metals, the presence of a superconducting 
contact adjacent to GMR multilayer does not suppress the GMR amplitude in the CPP 
geometry, except for special values of the parameters of the system. We think that such a 
particular situation has been considered in Ref. |J. The value of GMR depends not only 
on the spin-polarization of the electrons in the ferromagnetic layer, but also on the band 
parameters in the superconductor (in our case on value of R 2 ). Of course, since we used 
a simplified model, our results have only a qualitive nature. A more detailed analyzis of 
the dependence of the GMR on the microscopic parameters (exchange splitting, difference 
in T and j spin mean free paths), taking into account F/S interfacial scattering, will be 
presented in a forthcoming paper. 
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